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' Hypergeometric functions provide a useful representation of Feynman diagrams occuring in precision phenomenology. 

In dimensional regularization, the e-expansion of these functions about d = 4 is required. We discuss the current status 
+^ ' of differential reduction algorithms. As an illustration, we consider the construction of the all-order e-expansion of 

Othe Appell hypergeometric function F± around integer values of the parameters and present an explicit evaluation of 
the first few terms. 

\o : 

It is well-known that the hypergeometric representation is one of the most fruitful tools in the investigation of 
' analytical properties of Feynman diagrams, and for efficiently evaluating these diagrams. The manipulation of 
hypergeometric functions can be separated into two distinct problems: 

(A) deriving relations between hypergeometric functions with different parameters; 
■ (B) finding an algorithm for the construction of the Laurent expansion of hypergeometric functions. 

The first problem was solved by mathematicians [lj], wheres the second one, which has been analyzed mainly by 
physicists, is still a subject of active research. The series of our recent papers has been devoted to the latter problem. 
HH 3 There are three different ways to describe special functions of the type occurring in Feynman diagrams: 
(i) as an integral of the Euler or Mellin-Barnes type; 
CN| ' (ii) by a series whose coefficients satisfy certain recurrence relations; 

(iii) as a solution of a system of differential and/or difference equations (holonomic approach). 
qq For functions of a single variable, all of these representations are equivalent, but some properties of the function may 
be more evident in one representation than another. These three different representations have led physicists to three 
^ , separate approaches to developing the e-expansion of hypergeometric functions in Feynman diagram calculations. 

The most impressive result in the Euler integral representation was the construction of the all-order e-expansion of 
Gauss hypergeometric functions with special values of parameters in terms of Nielsen polylogarithms Q . Such Gauss 
hypergeometric functions are related to one-loop propagator-type diagrams with arbitrary masses and momenta, 
two-loop bubble diagrams with arbitrary masses, and one-loop massless vertex-type diagrams. 

The series representation has also been very useful and intensively studied. The first results of this type were derived 
in context of the so-called "single scale" diagrams Q . Particularly impressive results involving series representations 
were derived in the framework of "nested sums" in Ref. Q and in a generating function approach in Refs. 0, 

The differential equations satisfied by hypergeometric functions provide another approach. Results for Gauss 
hypergeometric functions expanded about integer and half-integer values of parameters were presented in Ref. 0, > 
and results for a special type of rational coefficients in Ref. 0] , while results for generalized hypergeometric functions 
p-Fp-i with integer values of parameters were presented in Ref. [4|. The iterated solution approach provides an 
efficient algorithm for constructing the e-expansion, since each new term is related to previously known terms. 

An important tool for the iterated solution is the iterated integral defined for each k by Iiflk, ctfc-i, ■ • ■ , flti; z) = 
Jo t-a k -^( a fc-ii ■ ■ ■ : a ii t)- A special case of this integral, 

Gm k ,m k - 1 ,...,m 1 (ak, ...,ai;z) = I( 0, . . . , ,ak,'"> 0,...,0 , av,z) , 
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is related to multiple polylogarithms 10 1 
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where we have used I(ai, . . . , a^; z) — I , 

As an illustration of our technique, we shall consider the Appell hypergeometric function F\ defined by the series 
F\(a, bi, hi , c; 21,22) — Em=oE^=o (shi+^lh^lhi For this function, the differential reduction algorithm |l| 

is equivalent to the statement that any function F±(a, b\, b 2 , c; z\, 22) can be expressed as a linear combination of 
function with arguments that differ from the original ones by an integer, Fi(a+mo; 61+mi; b 2 +m 2 ; c-Hm.3; z\, z 2 ), and its 
first derivatives R^F\{A+m; z\, 22) = J2j=i 2 Rj@j + Ro Fi(A; Z\, Z2), where 6 r = z r -^- ,r = 1,2, A = (a,bi,b 2 ,c) 
is a list of parameters, the m are lists of integers, and the Ri are polynomials in parameters a,b r ,c and z r . of basis 
functions The basis function ojq(zi, Z2) = Fi(ae, b\E, b 2 e, 1 + ce; Zi, z 2 ) satisfies the system of differential equations 
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Due to the analyticity of F\ with respect to its parameters, Eqs. (|T|) — (J2I> hold for every coefficient function coq (21 , 22) 
in the expansion ^0(21, 22) = 1 + Y]TL-t u>q (z\, z 2 )e k , and the coefficient equations are 
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Z r —u$\z 1 ,Z2) = uj ( r J \zi,z 2 ) ,r = 1,2 , 



uj{ j - l) (z 1 ,z 2 ) + b 1 uji j - 1) {z 1 ,Z2)+ab 1 ujV- 2 \z 1 ,Z2) , (3) 
4 j ~ 1) (z 1 ,Z2) + b 2 uj[ j - 1) (z 1 ,Z2)+ab2^- 2 \z 1 ,Z2) , (4) 

(5) 



where we have introduce new functions u> r (zi,z 2 ) — r u>o(zi, z 2 ), r — 1,2, which have e-expansions of the form 
uj r (zi,z 2 ) = 'Y^ = i^r k \zi 1 z 2 )e k , r = 1,2. To solve this system of first order differential equations, the boundary 
condition must be specified. Our choice is 

ctb 

ui r (z, 2) = z- — - — 2 Fi(l + ae, l + (&i + b 2 )e; 2 + ce; 2) and uj (z, 2) = 2 Fi(ae, (bi + b 2 )e; 1+ce; 2), 
1 + ce 

where the all-order e-expansion of the Gauss hypergeometric function has the form 0, Q 

2 F 1 (1+Ae,l + Be;2+C£;z)= ^ s ' c (2)e J ' 2F\(Ae, Be; 1 + Ce; z) = l + ABe 2 ^ W^ B ' C (z)e j , 



where p^' B ' C (2) and W°" Dl+ ° 2 ' c (z) are expressible in terms of generalized polylogarithms. [ldj j In particular, for the 
coefficient functions w?' (z±, z 2 ), we have 

4 j \z,z)+4 j \z,z) = a(b 1+ b 2 )p^ 2 +b ^ c (z) , 4 j \z,z) = a(b 1+ b 2 )W;i b 2 ^- c (z) . 
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The solution of Eqs. ©>© can now be written as 
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which should be supplemented by a self-consistency condition, 

[ f^- 1 \t,z 2 ) + h^- 1 \z 1 ,t)\ =a(b 1 + b 2 ) 
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where ^ (21,22) and ^(^i,^) are symbolic notation for the r.h.s. of Eqs. ([3]), respectively. The last 



equation ([5]) may be rewritten as (zi jj^ + z 2-g^ 
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(zi,z 2 ). The solution may be 



obtained by the method of characteristics, and we obtain a first-order ordinary differential equation depending on 
parameter C 



?' ^l^T^o" — u i ( z i, Czij+oj^' [z\, Czi). Collecting all of these relations together provides an 
iterative algorithm for constructing the all-order e-expansion of Appell hypergeometric function F\ around integer 
values of the parameters. Due to the fact that all integral solutions contain only the factors and j, we could expect 
that the result of iteration will always be expressible in terms of a particular case of multiple polylogarithms. [10( 
As an illustration, let us explicitly evaluate the first few coefficients. The first non-trivial term of the e-expansion 



corresponds to j = 2(u)r = 0): luq Z, (zi, z 2 ) = J2r=i 2 a ^r ^2 (z r ) and Ur"'(z\, z 2 ) = — J2r=i 2 a ^r l n (l — 2 r ) 
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In summary, the hypergeometric function approach to Feynman diagrams provides a promising route to organizing 
and calculating the higher-order processes which are becoming increasingly important in precision phenomenology. 
Conversely, mathematicians may find that results motivated by high-energy physics will uncover new relationships 
among these classes of functions, providing a fertile area for interaction between mathematics and physics. 
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